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Abstract. Bethe ansatz formulation is presented for several explicit examples of quasi 
• exactly solvable difference equations of one degree of freedom which are introduced recently 

by one of the present authors. These equations are deformation of the well-known exactly 
. solvable difference equations of the Meixner-Pollaczek, continuous Hahn, continuous dual 

Hahn, Wilson and Askey- Wilson polynomials. Up to an overall factor of the so-called 
pseudo ground state wavcfunction, the eigenfunctions within the exactly solvable subspace 
are given by polynomials whose roots arc solutions of the associated Bethe ansatz equations. 
The corresponding eigenvalues arc expressed in terms of these roots. 
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S '. 1 Introduction 

OO 

O 

Bethe ansatz method is one of the well-known solution methods for exactly solvable quantum 
systems (see, e.g. [T]) as well as for various spin models and statistical lattice models. In recent 
years the concept of exactly solvable quantum systems was drastically enlarged to include various 
examples of the so-called 'discrete' quantum mechanical systems [21 El SI O [T] , in which the 
Schrodinger equation is a difference equation instead of differential. Known examples of exactly 
solvable 'discrete' quantum mechanics are deformations of exactly solvable quantum mechanics, 
in which the momentum operators appear in exponentiated forms instead of polynomials in 
ordinary quantum mechanics. Their eigenfunctions are (g-)Askey scheme of hypergeometric 
orthogonal polynomials [3(9], which are deformations of the classical orthogonal polynomials, 
e.g. the Hermite, Laguerre and Jacobi polynomials. It is interesting to note that these examples 
are exactly solvable both in the Schrodinger and Heisenberg pictures [H [5]. That is, for each 
Hamiltonian of these examples, the Heisenberg operator solution for a special coordinate called 
the ' sinusoidal 1 coordinate can be constructed, as well as the complete set of the eigenvalues and 
the corresponding eigenfunctions. The eigenfunctions consist of the above-mentioned (g-)Askey 
scheme of hypergeometric orthogonal polynomials in the 'sinusoidal' coordinate. See [6j [7] for 
a comprehensive introduction of the 'discrete' quantum mechanics and the recent developments. 



*This paper is a contribution to the Proceedings of the 5-th Microconference "Analytic and Algebraic Me- 



thods V". The full collection is available at http://www.emis.de/journals/SIGMA/Prague2009.html 



2 



R. Sasaki, W.-L. Yang and Y.-Z. Zhang 



The domain of the Hamiltonian, its hermiticity (self-adjointness), the exact Heisenberg ope- 
rator solutions, the creation-annihilation operators and the dynamical symmetry algebras are 
explained in some detail. 

In this paper we present Bethe ansatz formulations and solutions for a family of Quasi- Exactly 
Solvable (QES) difference equations, which were recently introduced by one of the present au- 
thors |10| 111]. One of the main purposes of the present paper is to provide a good list of explicit 
Bethe-ansatz equations for the quasi-exactly solvable difference equations (|2.1ip . (|3.1ip . (|3.15p . 
(|4.9p . (|4.17p . and (|5.5|) . The list will be helpful for future research in various contexts of math- 
ematical physics. A quantum mechanical system is called quasi-exactly solvable (QES), if only 
a finite number (> 2) of eigenvalues and corresponding eigenvectors can be obtained exactly 
[j~2l [13]. Among various characterisation/identification of QES systems [I2j [HI [TU [15], we 
promoted a simple view that a QES system is obtained by a certain deformation of an exactly 
solvable quantum system. As shown explicitly for ordinary quantum mechanics by Sasaki- 
Takasaki [15], the deformation procedure applies to systems of many degrees of freedom as well 
as for single degree of freedom systems. This is in sharp contrast to the sl(2,R) characterisa- 
tion [13], whose applicability is limited to essentially single degree of freedom systems. Recently 
the deformation procedure was applied to exactly solvable 'discrete' quantum systems of one 
and many degrees of freedom to obtain corresponding 'discrete' QES systems [101 [TTj 116] . In 
this paper we present Bethe ansatz solutions for these QES difference equations. 

There are only a very limited number of examples of (quasi) exactly solvable difference 
equations which have been solved so far by the Bethe ansatz method. To our knowledge, 
the Bethe ansatz solutions were known only for some QES difference equations in connection 
with U q (sl{2)) [17], and for exactly solvable difference equations [I8j [191 [20] of the elliptic 
Rusijsenaars-Sclmeider system. Here we apply the Bethe ansatz formulation to several explicit 
examples of QES difference equations |10| 111] as deformation of exactly solvable 'discrete' quan- 
tum mechanics 0E], which are difference analogues of the well-known quasi exactly solvable 
systems, the harmonic oscillator (with/without the centrifugal potential) deformed by a sextic 
potential and the 1/ sin 2 x potential deformed by a cos 2x potential. As will be shown explicitly 
in the main text, these Bethe ansatz equations (|2.1ip . (|3.1ip . (|3.15p . (|4.9p . (|4.17p . and (|5.5p can 
be considered as deformations of the equations (|2.16p , (|4.14p , (|5.9p determining the roots of the 
corresponding (g-)Askey scheme of hyper geometric orthogonal polynomials [9] (the Meixner- 
Pollaczek, continuous Halm, continuous dual Hahn, Wilson and Askey- Wilson polynomials and 
their restrictions) which constitute the eigenfunctions of the undeformed exactly solvable quan- 
tum systems. The general structure of the quasi exactly solvable difference equations to be 
discussed in this paper and their solutions were explained in some detail in [101 ITT1 [To] , inclu- 
ding the domains of the Hamiltonians, the Hilbert spaces and hermiticity and the role played 
by the pseudo ground state wavefunctions. 

This paper is organised as follows. In Section [21 a QES discrete quantum mechanics is solved 
in the Bethe ansatz formalism. The Hamiltonian of the system is obtained by 'crossing' those 
of the Meixner-Pollaczek and the continuous Hahn polynomials as derived in [TT] . Section [3] 
provides the Bethe ansatz formulation of the QES discrete quantum mechanical systems which 
are deformations of the harmonic oscillator with a sextic potential as derived in [10] . The cor- 
responding eigenfunctions are deformations of the Meixner-Pollaczek and the continuous Hahn 
polynomials. Section [4] gives the Bethe ansatz solutions of the QES discrete quantum mechanical 
systems which are deformations of the harmonic oscillator with a centrifugal barrier and a sextic 
potential as derived in [10]. The corresponding eigenfunctions are deformations of the continuous 
dual Hahn and the Wilson polynomials. Section [5] offers a Bethe ansatz solution to a difference 
equation analogue of a QES system with the 1/ sin 2 x potential deformed by a cos 2x potential as 
derived in [10]. The corresponding eigenfunctions are deformations of the Askey-Wilson poly- 
nomials and their various restrictions [HI [9]. The final section is for a summary and comments. 
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2 Difference equation of the Meixner— Pollaczek type 

In this section we will discuss Bethe ansatz solutions for the discrete quantum mechanics ob- 
tained by deforming that of the Meixner-Pollaczek polynomials in [11]. To be more precise, 
the corresponding discrete quantum mechanics is obtained by crossing those of the Meixner- 
Pollaczek and the continuous Halm polynomials, that is, with the quadratic potential function 
of the continuous Hahn polynomial multiplied by a constant phase factor e~ 1 ^ of the Meixner- 
Pollaczek type. It was shown in [11] that this system is quasi exactly solvable. The corresponding 



Hamiltonian is: 

H d = V v ( x ) e ~ i9x V V ( X Y + V v ( x )* e+i9x V v ( x ) ~ v ( x ) ~ V ( X T + a M x (2.1) 

= A*A + a M x, a M d = -2Msin/3, M € Z+, (2.2) 

^ d = v^e-i 91 - A = e-i^y / V{xy-e^y / V{x), (2.3) 



V(x) = (ai + ix)(a 2 + ix)e lf3 , ai,a 2 eC, Re(ai) > 0, Re(a 2 ) > 0. 

It should be noted that the Hamiltonian is no longer positive semi-definite but the hermiticity 
is preserved. Let us introduce the so-called pseudo ground state wavefunction <j)o(x) \10\ 111]: 

fo{x) d = e^JTia! + ix)T(a 2 + ix)r(a* - ix)T{a* 2 - ix), 

as the zero mode of the A operator (|2.3p . A(fio = 0. The similarity transformed Hamiltonian TC 
in terms of t/>o, 

U = (j)" 1 oHo(j)v = V(x) (e"^ - 1) + V{x)* (e id * - 1) + a M x (2.4) 
= (ai + ix)(a 2 + ix)e- i(3 (e~ ldx - l) + (aj - ix){a* 2 - ix)e i(3 (e id * - l) - 2A^sin/3rr, 

acts on the polynomial part of the wavefunction. In the exactly solvable Meixner-Pollaczek case 
with V(x) = (a + ix)e~ l @ and the continuous Hahn case with V(x) = {a\ + ix)(a 2 + ix), the 
eigenfunctions are of the form 

(f) (x)P(rj(x)) 

in which P(r](x)) is a polynomial in 

rj(x) = x. 

After the deformation, it is obvious that 7i maps a polynomial in rj(x) = x into another and it 
is easy to verify 

Tix n = 2{—M. + n) sin/?x n+1 + lower order terms, n € Z + . 

This means that the system is not exactly solvable without the compensation term, but it is 
quasi exactly solvable, since TL has an invariant polynomial subspace of degree M: 

HV M C Vm, (2.5) 
V M = f Span [l,x, x 2 ,..., x M ] , dimV^ = M + 1. (2.6) 

Let ^(x) be one of the eigenfunctions of TL and E be the corresponding eigenvalue: 
WB(x) = EV(x), 
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namely, 

(a* - ix)(a* 2 - ix)e if3 (^(x + i) - <£(x)) 

+ (01 + ix)(a 2 + ix)e~ i/3 (V(x - i) - *(x)) - 2M sin (5x^(x) = EV(x). (2.7) 

Equations (|2.5p and (|2.6p imply that the eigenfunctions in the subspace Vm have the following 
form 

M M 

V(x) = J] (r?(x) - r)(xi)) = JJ(x - x,), (2.8) 

Z=l (=1 

where {x^ \ l = 1, . . . , A4} are some parameters which will be specified later by the associated 
Bethe ansatz equations (|2.1ip below. Substituting the above equation into (|2.7p and dividing 
both sides by ^(x), we have 

M _ _ ■ M _ _L ■ 

£ = -V(x) -V(xY -2Msm[3x + V{x)T\?—^ — - + V(xY TT £__?L±1. (2.9) 

- LJ - a; — X; - LJ - X — X; 

The r.h.s. of (|2,9p is a meromorphic function of x, whereas the l.h.s. is a constant. To make 
them equal, we must null the residues of the r.h.s. It is easy to see that the singularities of the 
r.h.s. only appear at x = Xj, j = 1, . . . , M. and x = oo. The residues at x = Xj vanish if the 
parameters {xj} satisfy the following Bethe ansatz equations 

TT xj -xi-i = VjxjT yjxj + j) -jjjxj) 
fj- Xj - xi + i V(xj) 7](xj) - r](xj - i) 

= ^^ W-«x i )(a S -ix i ) j = li _ }Mm (2 . u) 
V (Xj) (ai + iXj)(a 2 + ix j) 

Throughout this paper we use the complex conjugate potential function V(x)* in the 'analytical' 
sense in x, that is, for complex x 

V(x) = (ai +ix){a 2 + ix)e~ il3 , V(x)* = (a\ - ix){a* 2 - ix)e if3 , x e C. 

This convention is necessary for the above two equations (|2.1U|) . (|2.1ip to be valid, since the 
Bethe roots {xj} are in general complex. One can check that the r.h.s. of (|2.9p is indeed regular 
at x = oo, i.e., the residue at x = oo vanishes. By the Liouville theorem the r.h.s. of (|2.9|) 
is a constant provided that (|2.1ip is satisfied. One can get the value of the corresponding 
eigenvalue E by taking the limit of x — > oo for the r.h.s. of (|2.9p . Here we present the result: 



M 

E = M{M — 1) cos/3 + M((ai + a 2 )e" i/3 + (aj + a* 2 )e ip ) + 2 sin (3 V x h (2.12) 



where {x/} satisfy the Bethe ansatz equation (|2.1ip . The final term can be written as 
M 

2sin/3^??(x/). (2.13) 
l=i 

The wavefunction ^(x) (|2.8p becomes the eigenfunction of 7i in the subspace Vm (|2.6p provided 
that the roots of the polynomial ^(x) (|2.8p are the solutions of (|2.1ip . and then the corresponding 
eigenvalue is given by (|2.12|) . Since all the roots {x;} are on the same footing, it is natural that 
the eigenvalue E depends on the symmetric combination of them (|2.13p . 
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A few corollaries ensue from these results. For the special case of (3 = 0, the Hamilto- 
nian (|2.1|) is exactly solvable and the corresponding eigenvectors are related to the continuous 
Hahn polynomials. In fact, we obtain from (|2.12p 

\imE = M(M +ai + a 2 + a? + at - 1), (2.14) 

which is the eigenvalue corresponding to the degree M. continuous Hahn polynomial [2j [3l [6j. 
The Bethe ansatz equation ([2. lip now determines the zeros of the continuous Hahn polynomial. 
For a\, a 2 € K + , and a 2 — > oo limit, the above Hamiltonian H (|2.ip divided by a 2 reduces to an 
exactly solvable one corresponding to the Meixner-Pollaczek polynomials [2j [3j [6J . Correspon- 
dingly the eigenvalue formula (|2.12p gives 

lim E/a 2 = 2Mcos/3. (2.15) 

122^00 

The latter is the eigenvalue of the degree A4 Meixner-Pollaczek polynomial [6j. Note that the 
parameter (3 is related to the standard parameter <j) of the Meixner-Pollaczek polynomials as 
P = l£ — 4>. The corresponding Bethe ansatz equation 

M . / ■ \ 

Xj — Xl — l _ ^113 \ a l lX 3 



II j ~ Xl ~ l . =e^p—^l, J = 1,...,M, (2.16) 
ff-. Xj - xi + i \a\ + ix j) 

1+3 



now determines the zeros of the degree Ai Meixner-Pollaczek polynomial. As shown in § 4 of [2] , 
this equation reduces to that determines the zeros of the Hermite polynomial in an appropriate 
limit. 

3 Difference equation analogue of harmonic oscillator 
deformed by sextic potential 

There are two types of difference equations which are difference analogues of the sextic potential 
Hamiltonian [10]. The Hamiltonian is given by 

H = y/V(x)e- id *^V(xy + ^V{xr^ dx ^W)- (V(x) +V(x)*)+a M (x), (3.1) 

Type I: V(x) d = (a + ix)(b + ix)V (x), V Q (x) = c + ix, a,6,cGM+, (3.2) 

OiM{x) = f 2Mx 2 , 
Type II : V(x) d = (a + ix)(b + ix)V (x), 

Vq(x) = (c + ix)(d + ix), a,b,c,d 6l + , (3-3) 

a M (x) d = M (M - 1 + 2(a + b + c + d)) x 2 . (3.4) 

Here as usual V{x)* is the 'analytical' complex conjugate of V(x). 



3.1 Type I theory 

Here we will consider the difference equation of type I, while the difference equation of type II will 
be given in the next subsection. If V is replaced by Vq in (|3.2p and the last term in (|3.ip , (x) , 
is removed, 7i becomes the exactly solvable Hamiltonian of a difference analogue of the har- 
monic oscillator, or the deformed harmonic oscillator in 'discrete' quantum mechanics OEj. Its 
eigenfunctions consist of the Meixner-Pollaczek polynomials, with a special phase angle (5 = 0, 
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which is a deformation of the Hermite polynomials [21 [31 [21] . The quadratic polynomial factor 
(a + ix)(b + ix) can be considered as multiplicative deformation, although the parameters a, b 
and c are on the equal footing. On the other hand one can consider it as a multiplicative 
deformation by a linear polynomial in x: 

V(x) = (a + ix)Voi(x), Vbi(x) = f (b + ix)(c + ix), 

with Voi describing another difference version of an exactly solvable analogue of the harmonic 
oscillator [21 [3]. Its eigenfunctions consist of the continuous Hahn polynomials. 

Next let us introduce the similarity transformation in terms of the pseudo ground state 
wavefunction (j>o(x) as the zero mode of the A operator (|2.3p . A<fro = 0: 

4>q(x) = f y/T(a + ix)T(a — ix)T(b + ix)T(b — ix)F(c + ix)T(c — ix), 

H = ^ 1 oHofo = V{x)(e- i9 *-l)+V{x)*(e id *-l)+2Mx 2 . (3.5) 
Since the parity is conserved, that is 

H\x-+-x = H > 

it is easy to verify the action of the Hamiltonian Ti (|3.5p on monomials of x: 

( [n/2+l] 

a nd x n+2 ~ 2j , n<M-2, a n j G R, 

3=0 
[M/2] 

J2<jX M ~ 2j , n = M, a' n j G R. 

3=0 



Hx n 



(3.6) 



Here [m] is the standard Gauss' symbol denoting the greatest integer not exceeding or equal 
to m. According to the parity of the polynomials, there are two types of invariant subspaces Vm 
of H: 

UVm^Vm, (3.7) 
„ def f Sp&n[l,ri(x),...,r)(x) k ,...,r)(x) M / 2 ] , M : even, 

V -M = \ „ r , , , , T](x)=X , (3.8) 



dimV^vi 



x Span [l,r/(x), . . . ,rj(x) k , . . . ,r}(x)^ M 1 ^ 2 ] , M : odd, 

M/2 + 1, M : even, 
(A4 + l)/2, M:o&&. 



3.1.1 The case of even A4 

Let us introduce a positive integer TV such that M. = 2N. Equations (|3.7p and (|3.8p imply that 
the eigenfunctions of TL in the subspace Vm are °f t ne form 

JV Af 

^(») = II ^(x) - V(%i)) = Y[(x - xi)(x + xi). (3.9) 
i=i i=i 

Analogous calculation shows that the polynomial ^f(x) becomes the eigenfunction of Tt if the 
roots of the polynomial satisfy the Bethe ansatz equations 

tt {xj - xi - i){xj + xi - i) = Vjxj)* ^{xj+i) - rjjxj) 
f_}: (xj - xi +i)(xj + xi +i) V(xj) rj(xj) - rj(xj - i) 
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= (g ~ ix j)(P ~ ix j)( c ~ ixj)(2xj + Q j = \ N (3 11) 

(a + ixj)(b + ixj)(c + ixj)(2xj — i) ' 

The corresponding eigenvalue E is given by 

1 N 
E = -M(M - l)(M - 2) + (a + b + c)M{M - 1) + 2(a6 + ac + 6c)X - 4 V x 2 , (3.12) 

where {x;} satisfy the Bethe ansatz equations ([3. lip . Again the final term is symmetric in {xj} 
and can be written as —^Ylf=iV{ x j)- 

3.1.2 The case of odd M 

Let us introduce a positive integer N such that A4 = 2N + 1. (|3.7p and (|3.8p imply that the 
eigenfunctions of 7i in the subspace Vm are of the form 

N N 

= £ (r/(x) - T](xi)) = xY\_( x - xi)(x + xi). (3.13) 

Analogous calculation shows that the polynomial ^f(x) becomes the eigenfunction of Tt if the 
roots of the polynomial satisfy the Bethe ansatz equations 

(xj ~ tt ( x i ~ x i~ +xi-i) _ Vjxj)* (rjjxj +i)- rjjxj)) 
(xj +i) *£ (xj - xi+i)(xj +xi + i) V(xj) (r)(xj) - rj(xj - i)) 



(a — ixj)(b — ixj)(c — ixj)(2xj + i) 
(a + ixj)(b + ixj)(c + ixj)(2xj — i) 

The corresponding eigenvalue E is given by 



j = l,...,N. (3.15) 



1 - 
E = -M(M - l)(M - 2) + {a + b + c)M{M - 1) + 2(ab + ac + bc)M -A^xf, (3.16) 

i=i 

where {xi} satisfy the Bethe ansatz equations (|3.15l) . The expression for the eigenvalue is exactly 
the same as the even M. case (|3.12p . 

In this example, exactly solvable limits are also obtained by making one or two parameters 
go to infinity; for example a — > oo or both a — > oo and b — > oo. In the former case (a — > oo), 
the scaled Hamiltonian (|3.ip TL/a gives that of the continuous Hahn polynomials with real 
parameters b and c. The eigenvalue formulas (|3.12p and (|3.16p reduce to that of the continuous 
Hahn polynomials [21 [3J [6] : 

lim E/a = M(M + 2(b + c)-l). (3.17) 

a— >oo 

The Bethe ansatz equations (|3.1ip and ()3. 15j) are equivalent to (|2.1ip with (5 = and a\ = b 6 
EL)., a2 = c G R+. This assertion can be easily verified since the solutions of (|2.1ip are always 
paired {xj, —xj} including a zero Xj = for odd A4. In the latter case (a — > oo, b — > oo), the 
scaled Hamiltonian (|3.1|) Tt/(ab) gives that of the Meixner-Pollaczek polynomials with /? = 0, 
a = c. The eigenvalue formulas (|3.12p and (|3.16p reduce to that of the Meixner-Pollaczek 
polynomials [21 El [6] : 

lim E/(ab) = 2M. (3.18) 

a,b— >oo 

Again the Bethe ansatz equations (|3.1ip and ()3.15|) are equivalent to p. lip with (3 = and 
ai = c G R + . 
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3.2 Type II theory 

Another difference analogue of the sextic potential Hamiltonian has the same form as (|3.ip . 
with only the potential function V{x) and the compensation term a^(i) are different: 

V(x) d = (a + ix)(b + ix)Vo(x), Vq(x) = f (c + ix)(d + ix), a, b,c,d£ M.+ , 
a M (x) = M (M - 1 + 2{a + b + c + d)) x 2 . 

This Hamiltonian can be considered as a deformation by a quadratic polynomial factor (a + 
ix)(b + ix) of the exactly solvable 'discrete' quantum mechanics having the continuous Hahn 
polynomials as eigenfunctions [3j, another difference analogue of the harmonic oscillator. See 
the comments in Section 5 of |16j . 

The pseudo ground state wavefunction (j>o(x) is 

$0 (x) = f \/r(a + ix)T(a — ix)T(b + ix)T(b — ix)T(c + ix)T(c — ix)T(d + ix)T(d — ix). 

Again it has no node and it is square integrable. The similarity transformed Hamiltonian acting 
on the polynomial space is 

n d £ ^0710^ = V{x) (e~ id * - 1) + V{x)* (e id * - l) 

+ M{M -\ + 2{a + b + c + d))x 2 . (3.19) 



It is straightforward to verify the relationship (13. 6p and to establish the existence of the invari- 
ant polynomial subspaces. In the following subsections, we present Bethe ansatz solutions to 
eigenfunctions and eigenvalues of the Hamiltonian (|3.19p . 

3.2.1 The case of even Ai 

Let us introduce a positive integer N such that A4 = 2N and a polynomial function ^(x) of 
the form (|3.9D . Then ^(x) becomes the eigenfunction of TL (I3.19P with an even Ai: A4 = 2N 
provided that the roots of the polynomial satisfy the Bethe ansatz equations 

tt {xj-xi-i){xj+xi-i) = Vjxj)* (r](xj + i) - i]{xj)) 
jJL(xj-xi+i)(xj + xi + i) V(xj) (7](xj) -rj(xj -«)) 

(a — ixj)(b — ixj)(c — ixj)(d — ixj)(2xj + i) 
(a + ixj)(b + ixj)(c + ixj)(d + ixj)(2xj — i) ' 

The corresponding eigenvalue E is given by 

E = 2 E ( "V ) A ^ ~ (4A3 + {4M ~ 6)) E x *> ( 3 - 21 ) 

j=l V 3 J 1=1 

with the binomial coefficients 
'M\ Ml 



N. 



J J .i'iM ./)! 



0,...,A4, 



and the coefficients {Aj|j = 1, . . . ,4}, which are the elementary symmetric polynomials in the 
parameters {a, b, c, d} defined by 



4 

def 



V(x) = {a + ix)(b + ix)(c + ix){d + ix) = ^2Aj(ix) j , A 4 = 1. (3.22) 

j=0 
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3.2.2 The case of odd M 

Let us introduce a positive integer N such that A4 = 2N + 1 and a function ty(x) of the 
form (13. 13h . The polynomial ^f(x) becomes the eigenfunction of 7i (I3.19P with an odd A4: 
A4 = 2N + 1 provided that the roots of the polynomial satisfy the Bethe ansatz equations 

(xj ~ i) t~t (xj - xi - i){xj + Xj - i) _ Vjxj)* (rjjxj+i) - rjjxj)) 
(xj + i) fj- (xj - xi + i)(xj + xi + i) V(xj) (rj(xj) - rj(xj - «')) 

(a — ixj)(b — ixj)(c — ixj)(d — ixj)(2xj +i) . 
(a + ixj)(b + ixj)(c + ixj)(d + ixj)(2xj — i)' ^ 

The corresponding eigenvalue E is given by the same expressions as (|3.2ip - (|3.22p but with 
a different value of M: M = 2N + 1. 



4 Difference equation analogues of harmonic oscillator 
with centrifugal potential deformed by sextic potential 

The Bethe ansatz solutions for the difference equation analogues of the harmonic oscillator with 
the centrifugal potential deformed by a sextic potential [10] are discussed here. There are two 
types corresponding to the linear and quadratic polynomial deformations as discussed in |l(Jj . 
The corresponding exactly solvable difference equation has the Wilson polynomials [21 [3l [9] 
as the eigenfunctions. The Hamiltonians have the same form as (|3.ip . (|2.2p and (|2.3p . with only 
the potential function V(x) and the compensation term cim( x ) are different: 

Type I : V(x) = f (b + ix)V (x), a M (x) = f Mx 2 , 

Type II: V(x) d = (a + ix)(b + ix)V (x), (4.1) 

a M (x) = M (M - 1 + (a + b + c + d + e + /)) x 2 , (4.2) 
with a common Vq(x) 

Y/- / n def (c + ix){d + ix)(e + ix){f + ix) 

V (x) = ± 2ix(2ix + l) ' °' 6 ' C ' ^ + _ {1/2} - (4 ' 3) 

None of the parameters a, b, c, d, e or / should take the value 1/2, since it would cancel the 
denominator. Because of the centrifugal barrier, the dynamics is constrained to a half line; 
< x < oo. The type I case can also be considered as a quadratic polynomial deformation of 
the exactly solvable dynamics with Vq\{x): 

rr t if/ \ def , . w , . v T . Tr , , def {d + ix)(e + ix){f + ix) 

Typel: V(x ) = (b + ix c + ix )V 01 (x), V 01 (x ) = ± [\ \ 4.4 

2%x\lix + 1) 

which has the continuous dual Halm polynomials [2], El El Ej as eigenfunctions. This re-interpre- 
tation does not change the dynamics, since the Hamiltonian and A and A* operators depend 
on V{x). 

4.1 Type I theory 

Here we consider the difference equation of type I. The pseudo ground state wavefunction (po(x) 
is determined as the zero mode of the A operator (|2.3p . Acp = 0: 



Type I : <j) (x. 



J EI r(cy +ix)T( aj - ix) 
def Vj= 2 



y / F(2ix)T(-2ix) 
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in which the numbering of the parameters 

def def , def def , def def . / a r- \ 

a\ = a, a.2 = o, as = c, a<± = a, 05 = e, ag = j, (4.5) 

is used. It is obvious that (fio has no node in the half line < x < 00. 

The similarity transformed Hamiltonian acting on the polynomial space has the same form 
as before (E5 



TLx 



U = fa 1 oTLo(}) = V(x)(e- i9x - 1) + V{x)*(e idx - l) + Mr]{x). (4.6) 
Again the Hamiltonian is parity invariant, that is 

Although the potential V(x) has the harmful looking denominator l/{2ix(2ix+l)}, it is straight- 
forward to verify that TL maps a polynomial in rj{x) = x 2 into another, as 2ix+l oc r}{x—i)—r](x): 

n+l 

a n:j x 2n+2 - 2j , n<M-l, a n j € R, 

M 

J2<^ 2M ~ 2j , n = M, a' n>j ER. 
i j=o 

This is because Vq, which has the above denominator, keeps the polynomial subspace of any 
even degree invariant, reflecting the exact solvability. In other words, the exactly solvable 
discrete quantum mechanics corresponding to the undeformed potentials Vq{x) or Vbi(x) has 
the eigenfunction 

<f>o(x)P(r)(x)), r)(x) = x 2 , 

in which P{rf) is either the continuous dual Hahn polynomial (Voi(x), (|4.4p ). or the Wilson 
polynomial {Vq{x)^ (|4.3p ). This establishes that TL keeps the polynomial space Vm invariant, 

TLVm £ Vm> 

V m = f Span[l,x 2 ,...,x 2fc ,...,x 2A4 ], dimV^ = M + 1. (4.7) 

The above equations imply that the eigenfunctions of TL in the subspace Vm are °f the form 
M M 



1=1 1=1 

Analogous calculation shows that the polynomial ^f(x) becomes the eigenfunction of TL if the 
roots of the polynomial satisfy the Bethe ansatz equations 

TT {xj - xi - i)(xj + xi - i) = VjxjY y(xj+i) - rjjxj) 
(xj - xi + i)(xj + xi +i) V(xj) r](xj) - r](xj - i) 

= {b- ixj)(c - ixj)(d - ixj)(e - j Xj )(f - ixj) ^ j = i ,M (4.9) 
(b + ixj)(c + ixj)(d + ixj)(e + ixj)(f + ixj) ' 

Note that the kinematical factors ±2ix + 1 of V(x) and V(x)* are cancelled by rj(x ^fi) — r](x). 
The corresponding eigenvalue E is given by 



E = -M{M - \){M-2) + [b + c + d + e + f + -\ M{M - 1 
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M 

+ {b{c + d + e + /) + c{d + e + /) + d(e + /) + ef} M - J] (4.10) 

1=1 

where {xi} satisfy the Bethe ansatz equations (|4,9p . 

As expected, exactly solvable limits are obtained by two different ways; either one or two 
parameters go to infinity. In the former case, the scaled Hamiltonian (|3.1|) %/f gives that of 
the Wilson polynomials with four real parameters, b, c, d and e. The eigenvalue formula (|4.10p 
reduces to that of the Wilson polynomials [21 El E] 

lim E/f = M{M + b + c + d + e-l). (4.11) 

/-►oo 

In the latter case (e, / — > oo), the scaled Hamiltonian (|3.ip TL/{ef) gives that of the continuous 
dual Hahn polynomials with three parameters, b, c and d. The eigenvalue formula (|4.10p reduces 
to that of the continuous dual Hahn polynomials [2j El [6] 

lim E/(ef) = M. (4.12) 

e,/~>oo 

The Bethe ansatz equations (|4.9p in these limits determine the zeros of the Wilson and the 
continuous dual Hahn polynomials, respectively: 

f _> oo • TT (^j -i){ x j +xi-i) _ (b - ixj)(c - ixj)(d - ixj)(e - ixj) 
fJL (xj — xi + i)(xj + xi + i) {b + ixj)(c + ixj)(d + ixj){e + ixj)' 

e f -> oo • TT ( x j ~ g ~ + g ~ _ ( b ~ ^j)(c ~ ^j)(d - fgj) 
' (xj — x\ + i)(xj + + i) (6 + ixj)(c + ixj)(d + irr^) 

4.2 Type II theory 

Here we consider the difference equation of type II. The pseudo ground state wavefunction 4>o( x ) 
is determined again as the zero mode of the A operator (|2.3p . Afio = 0: 



Type II : 4>o(x) 



dcf V nj=i r ( a i + ix)T(aj - ix) 



v / T(2ix)T(-2ix) 



where the same numbering of the parameters as those in (|4.5|) has been used. The similarity 
transformed Hamiltonian acting on the polynomial space has the similar form as (|4,6p but with 
a different potential function V(x) (|4.ip and compensation term a^(i) (14.21) 



H = (pQ 1 oHo (f> = V(x) (e- id * - 1) + V(x)* {e id * - l) 

+ M(M-l + (a + b + c + d + e + f))x 2 . (4.15) 

It can be verified that the resulting Hamiltonian TC ()4.15j) keeps the polynomial space Vm 
defined in (|4.7p invariant. Like the type I case, this allows us to search the corresponding 
eigenfunction ^(x) of form ^>(x) = Yi{^i( x ~ x i)( x + x i)- Analogous calculation shows that the 
polynomial \&(x) becomes the eigenfunction of 7i (I4.15P provided that the roots of the polynomial 
satisfy the Bethe ansatz equations 

tt (xj - xi - i)(xj + xi - i) = VjxjY (rjjxj + i) - 7]( Xj )) 
i-fj ( X J ~ x i + i )( x j + x i + i ) V(xj) (r](xj) - rj(xj - i)) 
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= (g ~ ix j)( b - ixj)(c - ixj)(d - ixj)(e - ixj)(f - ixj) ^ M ^ ^ 

(a + ixj)(b + ixj)(c + ixj)(d + ixj)(e + ixj)(f + ' 

The corresponding eigenvalue £ is given by 

B = A 3 (*) + (2A 4 + A 5 ) (*) + 4(A 5 + l)(^ + 8(^)-(A s+! ( W -l))^. 

Here the coefficients {Aj|j = 1, . . . , 6} are the elementary symmetric polynomials in the para- 
meters {a, b, c, d, e, /} defined by 

6 

(a + ix)(b + ix)(c + + ix)(e + ix)(f + ix) = Aj(ix)- 7 , Ag = 1. 

3=0 

Various limits to the exactly solvable cases are almost the same as in the type I theory and will 
not be listed here. 

It is interesting to note that the type II Hamiltonian of Section [3l (|3.ip . (|3.3p . (|3.4p . is 
obtained from the type II Hamiltonian of Section HI (|3.1j) . (]4. If) . f)4.2|) as a formal limit e — > 0, 
/ -> V2: 

4Wsoction 4|e-»0,/-»l/2 = ^Section 3|tW-+27W- 

The Bethe ansatz equations together with the eigenvalue formulas are related in similar ways. 



5 Difference equation analogue of 1/ sin 2 a? potential 
deformed by cos 2x potential 

The last example is the difference analogue of the model discussed in Subsection 2.1.2 of |10j . 
1/ sin 2 x potential deformed by a cos 2x potential. In this case the corresponding exactly solvable 
difference equation has the Askey-Wilson polynomials [21 El El E] as eigenf unctions. The basic 
idea for showing quasi exact solvability is almost the same as shown above. 

As introduced and explored in [10] , this system is a quasi exactly solvable deformation of the 
exactly solvable dynamics which has the Askey-Wilson polynomials [21 [3j El [9] as eigenfunctions. 
Here we slightly change the notation from that of [10] for consistency with the rest of this paper. 
The range of the parameter x is now finite, to be chosen as 

< X < 7T, 

and we introduce a complex variable z and the sinusoidal coordinate rj(x): 
z = e lx , 7](x) d = cosx = (z + 2~ 1 )/2. 

def 

The unit of the shift is changed from 1 to a real constant 7 = logq, < q < 1. Then the shift 
operator e 7P can be written as 

e TP = e -*r& =q D D ^ 

dz 

whose action on a function of x can be expressed as z — > qz: 

e 1P f{x) = f(x - i 7 ) = q D f(z) = q z £f(z) = f(qz), with f(x) = f(z). 
Note that 7 < 0. 
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The Hamiltonian takes the form 



H = y/V(x) e^y/V(x)* + y/V(x)* e~"i p ^V{x) - V(x) - V(x)* + a M (x), 

= VV{x)q D ^V{xY + y/V(x)* q- ^/^ - (V(x) + V{xf) + a M {?) t (5.1) 
= A^A + a_A4(rc), aj^{x) = f —2abcdeq~ l (l — q M )i](x), 



At 



def 



-I 



(y/V(x)q% - VV(x)*q-^ , A d = i (q% ^V{x)* - q~% \/V(x)j , (5.2) 



, r ,,def n u , v dcf {I - bz)(l - cz)(l - dz)(l - ez 
V[x) = (1 - az)V {x), Vq{x) = 



(1 - z 2 )(l - qz 2 ) 

— 1 < a,b,c,d,e < 1. (5-3) 

The Hamiltonian is obtained by deforming the potential function Vq{z) by a linear polynomial 
in z. The parameter range (|5.3p could be enlarged to one real parameter (say, a) and two 
complex conjugate pairs (for example, b = c* , d = e*), but the absolute values must be less 
than 1, |o[ < 1, . . . , |e| < 1. 

The pseudo ground state wavefunction <j>o(x) is determined as the zero mode of the A opera- 
tor ([521), Acf) = 0: 



2.... 

oo 



(az,az 1 ,bz,bz 1 ,cz,cz 1 ,dz,dz 1 ,ez,ez~ 



1 ;q)oo' 



where (a±, . . . , a m ; 5)00 = f n^LiEKLo^ ~~ a jQ n )- Obviously <fro has no node or singularity in 
< x < 7r. We look for exact eigenvalues and eigenfunctions of the Hamiltonian (|5.ip in the 
form: 

H<t> = E(t>, 4>{x) = c/) (x)^(7](x)), rj{x) = cosx = (z + z _1 )/2, 

in which ^(rj(x)) is a polynomial in ij(x) or in [z + l/z)/2 = cosx. The similarity transformed 
Hamiltonian acting on the polynomial space has the form 



~dcf 



^0 1 W o (J) = V{z) (q D - 1) + V{z)* (q~ D - 1) - abcdeq" 1 (l - q M ) (z + M . 



Without the deformation factor 1 — az and the compensation term, the above Hamiltonian 7i is 
exactly solvable, that is, it keeps the polynomial subspace in rj(x) = [z + l/z)/2 of any degree 
invariant. The deformed Hamiltonian 7i is parity invariant Tt\ x ^_ x = Tt and it is straightforward 
to show the existence of an invariant polynomial subspace: 

V m = Span[l,»7(x),...,»7(x) fc ,...,»7(x) >( ] J dimV^ = M + 1. 

The above equations imply that the eigenfunctions of Tt in the subspace Vm are °f the form 
M M 



Vl/(x) = PJ(??( X ) ~~ v( x l)) = Y\( cosx ~ cosxi) = Pm{cosx). 



1=1 1=1 

Analogous calculation shows that *$>(x) becomes the eigenfunction of TC if the parameters {xi} 
satisfy the Bethe ansatz equations 

yr COsjXj - fry) - COS X[ = V(Xj)* TjjXj + jj) - T](Xj) 

~~ cos(xj + i'y) — cos X\ V(xj) r](xj) — r](xj — i-f) 
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(zj - a)( Zj - b){ Zj - c){ Zj - d){ Zj - e) j = l,...,M. (5.5) 



(1 — azj)(l — bzj)(l — czj)(l — dzj)(l — ezj)zj 

Note that Zj = e lXj and r](xj — z 7 ) = cos(xj — 27) = (qz + q~ 1 z~ 1 )/2, etc. Again the kinematical 
factors (1 — z ±2 )(l— qz ±2 ) of V(x) and V(x)* are cancelled by r](x^fij)—r](x). The corresponding 
eigenvalue E is given by 

E = (abed + abce + abde + acde + bcde)q~ 1 (q M - l) + q~ M - 1 

M 

- 2abcdeq M - 1 (l - q" 1 ) J^cosx,, (5.6) 

1=1 

where {xi} satisfy the Bethe ansatz equations (|5.5p . 

In this example, there are many ways to obtain exactly solvable dynamics; by making either 
one (e), two (d, e), three (c,d,e), four (b,c,d,e) or five (a,b,c,d,e) parameters vanish. The 
corresponding Hamiltonians are those describing the dynamics of the Askey- Wilson, continuous 
dual g-Hahn, Al-Salam-Chihara, big g-Hermite and q-Hermite polynomials, respectively [HI EJ 
[2j[3l[6]. The eigenvalue formula (|5.6p reduces to that of the Askey-Wilson polynomials for e = 

E = (q~ M - 1) (1 - abcdq M+1 ) , (5.7) 
and to a universal formula 

E = q~ M - 1, (5.8) 
for the rest [2 EJ [6] . The Bethe ansatz equations (|5.5j) for the restricted case of e = 



TT cos(xj -i 7 ) -cosx; = {zj - a){ Zj - b)( Zj - c)( Zj - d) ^ M 

j-J: cos(xj + ^7) — cosx/ (1 — azj)(l — bzj)(\ — czj)(l — dzj) ' 



then determine the zeros of the Askey-Wilson polynomials. Further restrictions d, e = 0, c, d, e = 
0, b, c, d, e = and a, b, c, d, e = determine the zeros of the continuous dual g-Hahn, Al-Salam- 
Chihara, big g-Hermite and g-Hermite polynomials, respectively [3 [6]. 



6 Summary and comments 

We have constructed Bethe ansatz solutions for the quasi exactly solvable difference equations of 
one degree of freedom introduced in [10] and [TT]. These quasi exactly solvable difference equa- 
tions are deformations of the well-known exactly solvable difference equations of the Meixner- 
Pollaczek, continuous Hahn, continuous dual Hahn, Wilson and Askey-Wilson polynomials. The 
eigenfunctions within the exactly solvable subspace are explicitly given by some polynomials 
(module a pseudo ground state wavefunction ^0) whose roots are solutions of the associated 
Bethe ansatz equations. The corresponding eigenvalues are expressed in terms of the solutions 
of the Bethe ansatz equations. These are difference equation counterparts of the results of 
Sasaki-Takasaki [15], which gave a Bethe ansatz formulation of the QES systems corresponding 
to the harmonic oscillator (with/without a centrifugal barrier) deformed by a sextic potential 
and the l/sin 2 x potential deformed by a cos2x potential. 

As in the exactly solvable quantum mechanics, the sinusoidal coordinates i](x) [H [5] play an 
essential role. The exactly solvable sector is spanned by it 



Span[l, i](x), t](x) k , rj(x) M ] , 
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then the Bethe ansatz equations take an almost universal form, ()2.10j) . (|3.10p . (|3.14p . (|3.20p . 
(|3.23|) . (|4.8|) . (|4.16|) and (|5.4|) . The limits or restrictions to the exactly solvable dynamics are 
demonstrated in detail including various eigenvalue formulas (|2.14p , (|2.15p , (|3.17j) , (|3.18p , (|4.1ip , 
(|4,12p . ()5.7p and (|5.8p . The Bethe ansatz equations reduce to those determining the zeros of the 
corresponding orthogonal polynomials, for example (|2.16p . (|4.13p . (|4.14p and (|5.9p . 

All the known quasi exactly solvable dynamics have the exactly solvable subspace consisting 
of polynomials in a certain variable (see, e.g. [22]). Our emphasis here is that the variable 
is the sinusoidal coordinate which plays the central role in the corresponding exactly solvable 
limits 

It should be mentioned that there exist some examples of deriving (quasi) exactly solvable 
difference equations in terms of Lie algebraic deformations of exactly solvable dynamics intro- 
duced by Turbiner and his collaborators [23]. These difference equations have shifts in the real 
direction ip(x±l) and the corresponding eigenfunctions have discrete orthogonality measures, in 
contrast to those discussed in this paper which have pure imaginary shifts, ip{x±i) or ip(x±i^), 
7 G R, and the corresponding eigenfunctions have continuous orthogonality measures. 

The deformations of exactly solvable dynamics for obtaining the quasi exactly solvable quan- 
tum systems introduced in [151 HUJ [IT] and discussed in this paper in detail, are of the simplest 
type, in which the compensation term is linear in the sinusoidal coordinate. Possibility of further 
deformations including quadratic compensation terms will be discussed in [24], in particular, for 
those quantum systems having discrete orthogonality measures [7]. 
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